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Abstract 

Effects of a flavor singlet vector interaction in the Polyakov - Nambu - Jona-Lasinio 
(PNJL) model are studied in combination with the axial U(l) breaking Kobayashi - 
Maskawa - 't Hooft interaction. Using a consistent cutoff scheme we investigate the 
QCD phase diagram and its dependence on the vector coupling strength g^. We find 
that the first order chiral phase transition at moderate baryon chemical potentials and 
its critical point, a generic feature of most NJL-type models without vector coupling, 
disappear for sufficiently large values of that are consistent with lattice QCD results 
at imaginary chemical potential. The influence of non-zero g^ on the curvature of the 
crossover boundary in the T — fi plane close to /i = is also examined. 



1. Introduction 

The equation of state and the phase diagram of dense, strongly interacting matter, 
both at high and low temperatures, are in the focus of much recent interest, from heavy- 
ion collisions to the physics of compact stars. Lattice QCD computations continue to 
struggle with the treatment of (real) baryon chemical potentials /i when their ratios with 
respect to temperature T exceed /i/T ~ 1. It is for this reason that investigations of the 
phase diagram at finite have so far relied mostly on models based on the symmetries 
and symmetry breaking patterns of low-energy QCD. 

Chiral models such as the Polyakov-loop-extended Nambu - Jona-Lasinio (PNJL) 
model [1, 2, 3, 4, 5] or the Polyakov- loop-improved quark-meson (PQM) model 0,0] and 
most recently, Dyson-Schwinger equation approaches Q , are frequently used to calculate 
(T, /i) diagrams at large fi/T. In the region of interest, such models usually display a first 
order chiral phase transition (i.e. a pronounced discontinuity in the chiral condensate, 
{qq){li,T)). More precisely, a chiral crossover starting from /i = at the transition 
temperature Tc ~ 0.2 GeV proceeds at /i > until it reaches a critical fic at which the 
transition becomes first order. Starting from this critical point, a first-order transition 
line extends downward in temperature and meets the /i axis (T = 0) typically at quark 
chemical potentials around /i ^ 0.3 GeV, or baryon chemical potentials /is — 3fi ^ 0.9 
GeV. The associated mixed phase covers a broad range of baryon densities ps from 
below pq — 0.16 fm""^, the equilibrium density of nuclear matter, to about three times 
Pq. However, this is the nuclear physics domain where matter is known to be composed of 
interacting nucleons rather than PNJL type constituent quarks, and no chiral first-order 
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phase transition is anywhere close. Hence the question must seriously be addressed Q 
about the existence of a first-order phase transition (other than the nuclear liguid-gas 
transition) at low temperatures and /is < 1 GeV, and consequently, about the existence 
of a critical point in the QCD phase diagram. 

The appearance of a first-order chiral phase transition in PNJL models is a character- 
istic feature of the simplest version of these models with a "classic" NJL type interaction 
between quarks [lol |. i.e. a chiral combination of scalar and pseudoscalar interactions. 
It has already been noticed in previous work 11, 12, [l^, [3, [Hj that the existence 
and location of the critical point in the phase diagram is extremely sensitive to addi- 
tional interactions in the effective Lagrangian, especially to vector couplings and terms 
representing the axial U{1) anomaly. In the present work, using a PNJL model, we 
systematically study the role of two interactions that are each capable of removing the 
first order transition completely from the phase diagram: a repulsive fiavor singlet vector 
interaction that shifts the baryon chemical potential, and the U(1)a breaking Kobayashi 
- Maskawa - 't Hooft interaction, the coupling strength of which also has a pronounced 
impact on the phase diagram. We focus in particular on the interplay of these two effects. 



2. The Model 



The present investigation is based on the three-flavor PNJL model llj 
mean field level. The NJL part of the Lagrangian for three quark flavors [1 
has the general form: 




dot 



(1) 



with the quark fields ip — {ipu, ^di^s) , the current quark matrix mg — diag {niu, md, mg) 
and the GcU-Mann matrices A" in SU{3) flavor space. We work in the isospin-symmetric 
limit with to„ = m^. For vanishing quark masses, the interaction terms in Eq.([T]), apart 
from Cdet, are invariant under the chiral J7(3)l x C/(3)i^ group. The U{1)a symmetry 
is broken by the axial anomaly in QCD. At the quark level, this symmetry breaking is 
introduced by adding the Kobayashi - Maskawa - t Hooft (KMT) interaction [2ll : 



Cdct = K [detf (1 + 75) V) + detf (1^ (1 - 75) 1^))] 



(2) 



with the KMT coupling strength K. Here detf are 3x3 determinants in flavor space. The 
combination of scalar-pseudoscalar, vector and axial vector interactions arises naturally 
if one starts from a QCD-inspired color current-current interaction and then performs a 
Fierz transform into color-singlet channels. In this case the coupling strengths have fixed 
values, Gv = Ga = G/2, relative to that of the scalar-pseudoscalar coupling G. 

In the mean field approximation, dynamical quark masses are generated by the gap 
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equations 



K 

Ms = TO, " + ^'^^ ' 

with the chiral (hght quark) condensate au = G{^/juipu) = o-d and the strange quark 
condensate as = G{tps4's)- At non-zero quark densities Uq = (V'V), the flavor singlet 
term of the vector interaction, the one involving ipjoX^ip, develops a non-zero expectation 
value while all other components of the vector and axial vector interactions have vanishing 
mean fields. In the present context we therefore focus on the reduced NJL Lagrangian 



-dot 



a=0 

and study the effects of the flavor singlet vector interaction 

(5) 

with varying coupling strength = |Gy. Note that the presence of such an interaction 
with positive acts repulsively and shifts the quark chemical potential, fi ^ — g^Uq, so 
as to reduce the baryon density at given temperature T . We note at this point that a non- 
zero and positive vector coupling strength close to g^ ~ 0.4 G was found to be important 
in reproducing two-flavor lattice QCD data for the phase diagram at imaginary chemical 
potential within the frame of a non-local PN JL model 2^ 2J1 ■ 



The U{\)a breaking KMT interaction is responsible for the large mass of the 77' meson 
relative to the remainder of the pseudoscalar meson octet. Indications of a possible in- 
medium rj' mass reduction [26*1 would suggests that K might vary as a function of T and 
ji. We will treat X as a constant but study results for different fractions K/Kq < 1, with 
Kq being the vacuum value of the U{1)a breaking coupling strength. 

Amongst several similar NJL parameter sets available in the literature [l^, [3 : all 
adjusted to reproduce primarily physical masses and decay constants of the pseudoscalar 
meson octet, we use here the set from We have checked that results with other sets 
differ in general only marginally. The parameters determined by fitting meson properties 
in vacuum are TOu = toj — 3.6 MeV, TOs = 87 MeV (consistent with current quark 
masses reported in 5^ at a renormalization scale of 2 GcV), a three-momentum cutoff 
A = 750 MeV, GA^ w 3.6 and KqA^ « 8.9. The suggested value in fisj for g^/G is 0.6 - 
0.7, so that a meaningful window for variations of the vector coupling is 0.4 < g^/ G < 0.7. 

The implementation of Polyakov loop dynamics in the PNJL model used here proceeds 
as in [4]. Quarks move via minimal gauge coupling in a homogeneous background color 
gauge field. The Polyakov loop itself is a path-ordered Wilson line winding around the 
imaginary time direction. The expectation value $ of the Polyakov loop, the order 
parameter for the conflnement-deconfinement transition, is a real number, = at 
mean field level. Its behavior is determined by the Polyakov loop effective potential, 



w($,$*,r) = -i 



oo + ai 1 — I + a2 I — 
3 



<i>*$ + 64(^) log [J ($,$*)]. (6) 



with the Haar measure term 



J ($, 6$*$ + 4 (^$*^ + - 3 ($*$)^ 



(7) 



Parameters are determined by fits to pure gauge lattice data with Tq = 270 MeV [J|: 
ao = 3.51, fli = —2.47, 02 = 15.2, 64 = —1.75. This Z// is added to the fermionic thermo- 
dynamic potentiaL The thermodynamic potential density in mean field approximation 
(with a zero-point energy contribution Q,q subtracted) is then given by 
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The -Ejdp I; cr, /i, $) are the quasiparticle energy eigenvalues from the fermionic propa- 
gator matrix and U the Polyakov loop effective potential detailed above in The 
regularization of fl with a three-momentum cutoff A is motivated as follows. The first 
term on the r.h.s. is the thermal quark quasiparticle energy contribution from momentum 
regimes below the cutoff where the N JL interaction is "turned on" . There the quasipar- 
ticles propagate in a condensate background with dynamically generated effective quark 
masses given in (0]). For momenta above the cutoff A, the NJL interaction is "off', and 
the quarks propagate as free particles with their current quark masses under the sole 
influence of the Polyakov loop. In this high-momentum regime all couplings are set to 
zero in accordance with the basic premises of the NJL model, i. e. 



Ej ^ Ef'' = + ^2 ± ^ ; G = K = =0 for |p | > A (9) 

(with — /i for quarks, +11 for antiquarks). This procedure avoids unphysical behavior of 
the quark condensates and guarantees that the pressure reaches the Stefan-Boltzmann 
limit at high temperatures. Indeed, the chiral condensates vanish in the high temperature 
limit so that the quarks are free and carry just their current masses according to Eq.Q. 







Finally, the mean fields entering D, are determined by the self-consistent equations 

dndndn 

dau das 9$ 



= 



(10) 



and 



on 

-9^ = "^ = ^^^ 



(11) 



with the quark baryon number density, Uq, and the baryon density, ps- 



^ The nonlocal version of the PN JL model 0, |23| is a further improvement that avoids the momentum 
cutoff altogether by introducing a momentum-dependent dynamical quark mass M (p) that connects 
smoothly to the perturbative QCD limit as p — >■ 00. 
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3. Results 



It has already been known that a term hke (O lowers the position in temperature of 
the critical point with increasing Ov > 0, until it disappears completely from the phase 
diagram at a sufficiently large 11, 21\- In our case, this happens for > g"'* — 
0.38 G. Lowering K leads to a similar effect (T^]. With the parameter set used here, the 
critical point disappears for K < Kc — 0.71 Kq. 

We have studied the effects of variations of both parameters, g^ and K, on the T- 
fi and T-pB phase diagrams. First, consider the phase diagram in the T-/i plane for 
three choices of g^ in Fig. [TJ The case = displays the familiar picture of a chiral 
crossover terminating as a second order transition at a critical point. From this critical 
point downward a first order chiral phase transition extends to the fj, axis at T = 0. 
For = G/2 the first order phase transition has already disappeared and turned into a 
continuous crossover all the way down in temperature. For even larger vector coupling, 
(7v — G, the crossover pattern is shifted to higher chemical potentials. Note that what 
is represented here as a "crossover line" does not mean that chiral symmetry is restored 
in its Wigner-Weyl realization beyond that line. It just means that, along the line, the 
chiral condensate is reduced in magnitude to half of its value at zero temperature and 
H = 0. The crossover is in fact smooth and extends to densities more than four times 
the baryon density of normal nuclear matter. 
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Figure 1: PNJL phase diagrams for three differ- 
ent vector coupling strengths: (a) = 0, (b) 
Qv = G/2 and (c) = G. Dashed hues denote 
(Ju {T, /i) /cu (T = 0, = 0) = 0.5 at crossover 
transitions. The sohd hne is the first order transi- 
tion ending in the critical point. 



Figure 2: Curvature k of the crossover boundary 
near Tc as function of the vector coupling strength 
gv. Dots: results of calculations using the same in- 
put as in Fig. [T] Uncertainty measures attached 
to the dots give impression of the half width over 
which the crossover spreads. Upper horizontal 
line with narrow error band: (2-(-l)-flavor lattice 
QCD result of Ref . (28il : broad grey band: result of 
Ref.|2i|. 



Shifting the chiral crossover to regions of larger /i implies a flattening of the curvature 
of the crossover boundary in the neighborhood of fj, — 0. This is shown in Fig. [21 The 
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curvature of the crossover boundary is defined by 



or equivalently, k = —T^ 



Te(0) ; ' ^ d^i 



(12) 



Standard NJL or PNJL calculations without vector interaction usually produce a curva- 
ture parameter n that is too large in comparison with lattice results. A recent accurate 
lattice QCD computation [2^ gives n ~ 0.06. Evidently, a reasonably large vector cou- 
pling strength is capable of approaching such a small curvature. 




Figure 3: Phase structure in the PNJL model displayed as a function of the vector coupling strength 
(relative to the scalar-pseudoscalar coupling G) and the U{1)a breaking (KMT) coupling K (relative to 
its vacuum value Kq). The diagonal line separates regions in which a first order chiral phase transition 
takes place from those with smooth crossover. The grey area indicates the (unphysical) domain in which 
the mixed-phase region of a chiral first-order transition would overlap with the coexistence region of the 
liqu id-gas phase transition in nuclear matter as calculated using in-medium chiral effective field theory 



Next, consider the phase structure in scenarios with varying vector coupling and 
KMT coupling K . Increasing as well as reducing K has the effect of removing the 
first order chiral phase transition altogether from the phase diagram. The systematics 
is displayed in Fig. |31 Evidently, the appearance or disappearance of a first order phase 
transition is very sensitive to both the vector interaction and the strength of the axial 
C/(l) anomaly. For any value K < Kq, the 1st order transition turns into crossover when 
9v ^ 0.4 G. If K is reduced in the medium relative to the vacuum, the crossover scenario 
is initiated for even smaller values of gv/G. (We recall again that preferred values of 
gv range between 0.4 and 0.7 judging from lattice QCD thermodynamics at imaginary 
chemical potential and from earlier NJL phenomenology.) 

The first order chiral transition realized at — and viewed in the plane of tem- 
perature and baryon number density pe is known to run into an unphysical situation 
[9|: the phase mixing region associated with this 1st order chiral transition overlaps with 
the coexistence region of the liquid-gas phase transition in nuclear matter, realistically 
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calculated within the framework of in-medium chiral effective field theory [30| and prop- 
erly using nucleons rather than quarks as active fermionic degrees of freedom. Such an 
overlap, marked by the grey area in Fig. |31 is simply forbidden by well known facts of 
nuclear physics. This unphysical behavior is removed as soon as the repulsive vector 
interaction is introduced with reasonable values of its coupling strength, > 0.4 G, and 
the chiral transition turns to a smooth crossover. 



4. Summary and Discussion 

1. A first lesson from the present study is the strong sensitivity of the existence and 
location of a critical point to vector and U{1)a breaking interactions in any NJL or PNJL 
model. Realistic magnitudes of the vector coupling strength avoid the first-order chiral 
phase transition altogether and likewise dismiss the associated existence of a critical point 
in the phase diagram. 

2. While the exact value of the vector coupling strength is uncertain, it is likely 
to be larger than 0.4 times the strength of the chiral scalar-pseudoscalar interaction in 
NJL type models. This is sufficient to induce a smooth chiral crossover scenario at high 
densities and low temperatures, leaving well known nuclear physics constraints at lower 
baryon densities intact. 

3. A sufficiently strong flavor singlet vector interaction of the kind discussed in this 
work is likely to provide the necessary repulsion at high density that helps supporting two- 



solar-mass neutron stars [3l[ by making the equation of state of such systems sufficiently 
stiff m, 111]- Work along these lines is in progress. 

In summary, considering results from various sources using diverse methods, we tend 
to find a strong preference for a QCD phase diagram without a first order chiral phase 
transition, but rather with a smooth crossover at high densities and low temperature. 
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